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ABSTRACT 

We present recent contemporaneous X-ray and optical observations of the Be/X-ray binary system A 0535+26 with 
the Fermi /Gum ma-ray Burst Monitor (GBM) and several ground-based observatories. These new observations are 
put into the context of the rich historical data (since ~ 1978) and discussed in terms of the neutron-star-Be-disk 
interaction. The Be circumstellar disk was exceptionally large just before the 2009 December giant outburst, which 
may explain the origin of the unusual recent X-ray activity of this source. We found a peculiar evolution of the pulse 
profile during this giant outburst, with the two main components evolving in opposite ways with energy. A hard 
30-70 mHz X-ray quasi-periodic oscillation was detected with GBM during this 2009 December giant outburst. 
It becomes stronger with increasing energy and disappears at energies below 25 keV. In the long term a strong 
optical /X-ray correlation was found for this system, however in the medium term the Ha equivalent width and the 
V-band brightness showed an anti-correlation after ~2002 August. Each giant X-ray outburst occurred during a 
decline phase of the optical brightness, while the Ha showed a strong emission. In late 2010 and before the 2011 
February outburst, rapid V/R variations are observed in the strength of the two peaks of the Ha line. These had 
a period of —25 days and we suggest the presence of a global one-armed oscillation to explain this scenario. A 
general pattern might be inferred, where the disk becomes weaker and shows V/R variability beginning ~6 months 
following a giant outburst. 

Key words: accretion, accretion disks - pulsars: individual (A 0535+26) - stars: emission-line, Be - stars: 
neutron - X-rays: binaries 
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1. INTRODUCTION 

X-ray binaries are composed of a donor star, usually still 

on the main sequence, and an accreting compact object, which 
is either a neutron star or black hole. If periodic pulsations are 

detected from an X-ray binary, the compact object must be a neu- 

tron star and the system is called an accreting X-ray pulsar. The 

largest sub-group of high-mass X-ray binaries are the so-called 

Be/X-ray binaries (BeXRBs) in which the companion is a 

dwarf, subgiant, or giant OBe star. Be stars are rapidly rotat- 

ing objects with a quasi-Keplerian disk around their equator. 
The ultimate cause of the formation of the disk is still under de- 
bate (Townsend et al. 2004). In BeXRBs the optical and infrared 
emission is dominated by Be star companion and characterized 
by spectral lines in emission (particularly those of the Balmer 

series) and an IR excess, revealing the physical state of the mass 
donor component. In particular, the detection of emission-line 
features from the companion will confirm the presence of a disk 
around its equator (Coe et al. 2006). 

Most BeXRBs are transient systems. Historically, their out- 
bursts have been divided into two classes. In the long term, the 
X-ray variability of the transient BeXRB is characterized by 


type I (or normal) outbursts. These are regular and 
(quasi)periodic outbursts, normally peaking at or close to perias- 
tron passage of the neutron star, and reaching peak luminosities 
L x G 10 37 erg s . Type II (or giant) outbursts reach luminosi- 
ties of the order of the Eddington luminosity for a neutron star 
(i.e., when the gravitational attraction balances the outward ra- 
diation force on the accreting material; L x ~10 38 erg s l ; Frank 
et al. 2002) and become the brightest objects of the X-ray sky. 
During type II outbursts, an accretion disk may form. Unlike 
normal outbursts, giant outbursts have no consistently preferred 
orbital phase (Wilson et al. 2008). 

A 0535+26 is a transient BeXRB pulsar, discovered by Ariel V 
in 1975 (Rosenberg et al. 1975; Coe et al. 1975). The orbital 
period is ~ 1 1 1 days with a pulse period of ~ 103 s in a eccentric 
orbit (e ~ 0.47; Finger et al. 1996, and references therein). 
The optical counterpart is the 09.7 Hie star HDE 245770 
(Giangrande et al. 1980). An exhaustive review of observations 
at different wavelengths from 1970 until 1989 is found in 
Giovannelli & Graziati (1992). 

This system shows both giant and normal outbursts. The 
model of Okazaki & Negueruela (2001) attempts to explain 
the origin of type I and II outbursts by the truncation of 
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the circumstellar disk at a resonance radius between the disk 
Keplerian velocity and the orbit of the neutron star. Whether or 
not type I or II outbursts are produced depends on the eccentric- 
ity of the system. For the case of the moderate eccentricity of 
A 0535+26 (e = 0.47), both types of behavior are possible as 
the disk expands and contracts between different resonant radii 
(Haigh et al. 2004; Coe et al. 2006). We would therefore expect 
some correlation between the disk status as inferred from Ha 
equivalent width (EW) and profile and the presence or otherwise 
of X-ray outbursts. 

High-dispersion optical spectroscopic observations of this 
source during and after the 2009 December giant outburst 
showed drastic variabilities and indicated the existence of a 
warped component (Moritani et al. 201 1). A gas stream from a 
dense part of the Be disk to the neutron star is implied. Yan et al. 
(2012) found an anti-correlation between the optical brightness 
and the Ha intensity in 2009 observations. This indicated that 
a mass ejection event had taken place before the 2009 giant 
X-ray outburst. Near-IR monitoring during the 2011 February 
giant outburst and X-ray quiescent phases showed a ~ 1 2% 
reduction in the near-IR flux during the periastron passage of the 
neutron star (Naik et al. 2012). No changes in the spectra were 
found. 

During the 1994 giant outburst, broad quasi-periodic oscilla- 
tions (QPOs) from 27 to 72 mHz were detected (Finger et al. 
1996), confirming the presence of an accretion disk. In the 1989 
March/ April giant outburst, two cyclotron resonance scatter- 
ing features were detected at 45 keV and 100 keV (Kendziorra 
et al. 1994), from which a magnetic field of B ~ 4 x 10 12 G 
was inferred. First cyclotron line studies of the recent outbursts 
of A 0535+26 that took place in 2009/2010 have been pre- 
sented by Caballero et al. (2011b). The cyclotron line energy 
evolution revealed no significant variation of the cyclotron line 
energy. High-resolution grating spectroscopy of A 0535+262 
with Chandra detected for the first time an absorption consis- 
tent with the presence of a highly ionized outflow (Reynolds & 
Miller 2010), supporting the picture of a continued large-scale 
outflow at large accretion rates but in a non-jet form. Gamma-ray 
observations of A 0535+26 with VERITAS at very high energies 
(VHE; E ^ 100 GeV) were analyzed by Acciari et al. (2011) 
during the X-ray giant outburst in 2009 December. No VHE 
emission was evident at any time. They also examined data 
from the contemporaneous observations of A 0535+26 from 
the Fermi /Large Area Telescope at high-energy photons (E > 
0.1 GeV) and failed to detect the source at GeV energies. 

In this work we present contemporaneous X-ray/optical 
observations of the BeXRB system A 0535+26. We will interpret 
the data in terms of the B e-disk interaction with the neutron star 
companion. We include new optical observations from 2009 
October until 2011 March (MJD 55082 to 55665), from the 
Spanish Astronomical Observatories of Sierra Nevada (OSN), 
University of Valencia (OAO), and the Liverpool Telescope (LT, 
La Palma) as well as the Chinese Xinglong Station of National 
Astronomical Observatories (NAOC), and simultaneous X-ray 
observations by the Fermi/ GBM. 

2. OBSERVATIONS AND ANALYSIS 

In order to study the long-term X-ray behavior of A 0535+26, 
we have collected observations from a variety of sources. 
X-ray flux measurements by distinct missions from 1970 un- 
til 1989 were obtained from Giovannelli & Graziati (1992). 
Earth occultation observations with Compton Gamma Ray 


Observatory /BKTSYL were provided by the team 13 covering 
the period from 1991 to 2000. Rossi X-Ray Timing Explorer 
(RXTE) data provided by the All Sky Monitor (ASM) 14 team and 
Swift /BAT transient monitor results provided by the Swift /BAY 
team 13 extend the monitoring up to date. X-ray timing measure- 
ments were obtained from RXTE and GBM. Optical measure- 
ments were obtained from several sources, described below. 

In addition, in this work we present results for the continuous 
monitoring of A 0535+26 by Fermi GBM (Meegan et al. 
2009) since 2008 June 1 1 . The GBM is an all-sky instrument 
sensitive to X-rays and gamma rays with energies between 
~8 keV and ~40 MeV. GBM includes 12 sodium iodide 
(Nal) scintillation detectors and 2 bismuth germanate (BGO) 
scintillation detectors. The Nal detectors cover the lower part 
of the energy range, from 8 keV to about 1 MeV. The BGO 
detectors cover the energy range of ~ 1 50 keV to ~40 MeV. 
Only data from the Nal detectors were used in the analysis 
presented here. 

All the total-flux light curves for A 0535+26 by GBM were 
obtained using the Earth occultation technique 1 ' 1 on CSPEC data 
(128 energy channels every 4.096 s). Since GBM is not a pointed 
or imaging instrument, in order to determine fluxes for known 
sources, we measure the change in the count rate observed in the 
Nal detectors when the source enters or exits Earth occultation. 
A detailed description of the Earth occultation technique is given 
in Case et al. (2011) and Wilson-Hodge et al. (2012). Here we 
give a brief summary. For each occultation step, the count rates 
in a 4 minute window centered on the occultation time are fitted 
with a quadratic background plus source terms for the source of 
interest and each interfering source that appears in the window. 
The source terms consist of an energy-dependent atmospheric 
transmission model and a time-dependent model count rate, 
derived from the time-dependent detector response convolved 
with an assumed source spectrum. Each energy channel and 
each detector is fitted independently. For each source term a 
scaling factor is fitted, along with the quadratic background 
coefficients. When multiple detectors are included in the fit, the 
weighted mean for the scaling factor of the source of interest 
is computed for each energy channel. The mean scaling factor 
is then multiplied by the predicted flux in each energy channel 
to obtain flux measurements for the source of interest. In the 
case of A 0535+26, the assumed spectrum (Kendziorra et al. 
1994) consists of a power law with a high-energy cutoff and a 
Lorenzian cyclotron absorption model. Case et al. (2011) and 
Wilson-Hodge et al. (2012) find that flux values do not depend 
strongly on the choice of assumed model. Large deviations from 
expected fluxes are only seen if the model is very inconsistent, 
e.g., a soft model with a rapid turnover for a hard source or a 
very hard model for a soft source. Since each energy channel is 
fitted independently and channel-to-channel variations are not 
constrained by the model, hardness ratios (HRs) also appear 
not to strongly depend upon the choice of spectral model. 
A 0535+26 was not detected in the Nal detectors above 100 keV, 
even in the 10 days around the giant outburst peak. Since the 
Nal detector area is comparable to or larger than the BGOs 
(depending upon the combination of Nal detectors observing 
the source), the BGOs were not used. 

Timing analysis was carried out with GBM CTIME data (with 
eight channel spectra every 0.256 s) from channels 0, 1, 2, 

13 http://gammaray.msfc.nasa.gov/batse/occultation 

14 http://xte.mit.edu/ASM_lc.html 

15 http://heasarc.gsfc.nasa.gov/docs/swift/results/transients 

16 http://heastro.phys.lsu.edu/gbm 
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Figure 1 . Top: history of the X-ray outbursts undergone by A 0535+26 in arbitrary units. The peak intensities of the different outbursts are only intended to be 
illustrative of the times and types of events, since they have not been corrected for the different energy bands. Middle: long-term history of the evolution of the Ha 
equivalent width. Bottom: evolution of the visual magnitude. The compilation of Ha EW and visual magnitudes come from different authors described in the text (see 
Section 3.3). 

(A color version of this figure is available in the online journal.) 

and 3 (-8-12 keV, -12-25 keV, -25-50 keV, and 
~50-100 keV, respectively). After the appropriate Nal detector 
rates (from all 12 detectors and the four lowest channels) are 
selected, an empirical background model is fit and subtracted. 

The residuals to this, which contain the pulsations of the pulsars, 
are used to obtain pulse frequencies, pulse profiles, and pulsed 
fluxes. Then, a pulsed search was made over a narrow frequency 
band. A detailed explanation of our technique can be found in 
Finger et al. (2009a) and Camero-Arranz et al. (2010). 

The RXTE (Bradt et al. 1993) carries three instruments on 
board. The Proportional Counter Array (PCA; Jahoda et al. 

1996) is sensitive from 2 to 60 keV. The Fligh-Energy X-Ray 
Timing Experiment (Gruber et al. 1996) extends the X-ray 
sensitivity up to 200 keV. Monitoring the long-term behavior 
of some of the brightest X-ray sources, the ASM (Levine et al. 

1996) scans most of the sky every 1.5 hr at 2-10 keV. RXTE 
observed A 0535+26 since 2005. Recent observations since the 
beginning of 2009 have been analyzed in this work. For the 
pulse profile study and for each available observation, we have 
analyzed either GoodXenon or event-mode PCA data (time bin 
size 125 ms, 64M energy channels) using Ftools V6.9. The 
available PCA pulse profiles obtained in the ~8-14 keV and 


~ 14-21 keV energy bands are consistent with the GBM pulse 
profiles from channels 0 and 1 . 

Regarding the optical data, we include photometric archival 
observations from HIPPARCOS, ASAS, INTEGRAL /OMC, 
and recent ground-based observations from the 0.9 and 1.5 m 
telescopes at the OSN and the 60 cm telescope at the OAO. 
We transformed the HIPPARCOS filter to visual magnitude, 
following the procedure found in Harmanec (1998) and applying 
B—V and U — B values from Lyuty & Zaitseva (2000). We 
computed the V magnitude from Coe et al. (2006) based on 
the work done by Grundstrom et al. (2007). As we mentioned 
before, the spectroscopic observations come from the 2 m 
telescope at LT and the 2.16 m telescope at NAOC. Standard 
tools and photometry packages in IRAF were used for the 
reduction and analysis of the new data set of observations. 

3. OUTBURSTS HISTORY 
3.1. Long-term History 

Figure 1 shows a long-term X-ray/optical overview of 
A 0535+26/HDE 245770. The pre-BATSE X-ray flux mea- 
surements found in Giovannelli & Graziati (1992) come from a 
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Figure 2. Left: BAT light curve of the double-peaked normal outburst occurred in 2009 August (black circles). Overplotted is the spin frequency history during the 
event by GBM (red squares). Right: GBM 12-25 keV pulsed flux vs. time for the 2010 July normal outburst, with the frequency history also overplotted. 

(A color version of this figure is available in the online journal.) 


large variety of instruments, operating in different energy bands, 
and were presented only in Crab units. Due to the fact that the 
conversion of those measurements to real flux values is not 
straightforward and in order to show a global picture of this 
system, all the X-ray flux/count rates were converted into Crab 
units and then normalized using the maximum value for the 
whole data set (that is, BATSE). Therefore, the peak intensities 
of the different outbursts are intended to be illustrative of the 
times and types of outbursts. This will allow direct comparisons 
with the optical data and help us to understand the underlying 
Be-neutron-star interaction. 

The typical outburst peak luminosity for A 0535+26 during 
a normal outburst is ~400 mCrab and between 3 and 6 Crab 
for a giant one, in the 15-50 keV energy range. In the up- 
per panel of Figure 1, we see that after the giant outbursts in 
1975 (the discovery), 1977, 1980, 1983 (Sembay et al. 1990), 
and 1989, BATSE detected this source in 1994 (Bildsten et al. 
1997). The 1994 giant outburst was preceded by three normal 
outbursts and followed by two small normal outbursts. Then 
A 0535+26 went back to quiescence for almost 1 1 years. Only 
a few observations were reported during that period (Coe et al. 
2006; Hill et al. 2007; Negueruela et al. 2000; Orlandini et al. 
2004). The source renewed activity in 2005 May with another 
giant outburst and two normal outbursts. Each of those nor- 
mal outbursts began approximately a week before the perias- 
tron passage in the A 0535+26’s 111 day orbit. Unfortunately, 
the 2005 giant outburst was poorly observed due to the close- 
ness of the Sun to the source. There are only a handful of 
RXTE/ ASM dwells and serendipitous Swift /BAT observations 
near the outburst peak. In addition, a pointed RHESSI observa- 
tion (Smith et al. 2005) has been reported. 

3.2. Recent Outbursts 

After almost three years of being in quiescence A 0535+26 
was again active in 2008, undergoing four consecutive normal 
outbursts plus a giant outburst in 2009 December (see Figure 1). 
The very large event occurred in 2009 December was foreseen 
due to the fact that new X-ray activity was detected 19 days 
before periastron. Unexpectedly, another outburst peaking at ~ 1 
Crab (in the 15-50 keV band) occurred, with intermediate low 
X-ray activity being detected. This event was followed by two 
weak normal outbursts. In 201 1 February another giant outburst 
took place, approximately 1 1 days before the next periastron 
passage (2011 February 20). The outbursts of higher fluxes 
began at an earlier orbital phase. It is to be noted that this is the 


shortest interval in which two (or maybe three?) giant outbursts 
have been reported from A 0535+26. 

3.2.1. Some Peculiar Normal Outbursts 

Initial flaring of an outburst in 2005 September was seen 
as a double peak by RXTE. The highest flux in this outburst 
was reached during the short spike prior to periastron passage, 
and it was followed by a more gradual rise to a peak and then 
an exponential tail (Finger et al. 2006). Similar pre-periastron 
peaks have been seen in other BeXRBs such as 2S 1845-024 
and EXO 2030+375 (Camera Arranz et al. 2005; Wilson et al. 
2008). Using Swift /BPS data, Postnov et al. (2008) found that 
this was a regular normal outburst but exhibiting large flaring 
activity. A series of small flares taking place at the beginning 
of this event were particularly intense compared with the rest 
of the outburst. We also want to note that the remaining part of 
the outburst, as well as other outbursts of the source, presents a 
series of small flares instead of a smooth evolution of the flux 
with time. 

More recently, A 0535+26 underwent two double-peaked nor- 
mal outbursts in 2009 August and 2010 July (see Figure 2). We 
point out that the one in 2009 August showed two approximately 
equal peaks in both duration and intensity (see the left panel of 
Figure 2) without the typical shape of pre-periastron peaks de- 
scribed before. The last peculiar outburst in 2010 July was only 
entirely measured by GBM. One could note that this outburst 
is also noticeable in the BAT data, but these unfortunately have 
a gap up to the first peak (MJD 55374) because of the Sun (on 
2010 June 25, the Sun was ~8?5 from A 0535+26). Once more 
we found no resemblance with previous double-peaked out- 
bursts. The first peak was smaller, shorter, and totally detached 
from the main outburst peak (see the right panel of Figure 2). 
During the 2005 pre-outburst phase the pulse profile changed 
dramatically (Camero-Arranz 2007; Caballero et al. 2007), as 
did the value of the center of the cyclotron line energy (Caballero 
et al. 2008). However, in the 2009 and 2010 normal outbursts 
no changes have been found neither with RXTE / PCA nor GBM 
(see Section 5.3). 

3.2.2. Spectral Hardness 

Figure 3 shows the hardness-ratio-intensity (HRI) diagram 
analysis carried out for this source. This allows us to study the 
spectral variability of A 0535+26 without the supposition of 
any spectral model. GBM total-flux light curves from 10 keV 
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Figure 3. Hardness ratio analysis using GBM Earth occultation data for A 0535+26. The vertical axis indicates the flux in the 10-80 keV band. The HR has been 
defined by 20—40 keV/ 10-20 keV. The presence of observations with lower hardness ratio points to a change in the spectral continuum, i.e., a spectral cutoff at higher 
energies or an increased contribution from the softer components related to the neutron star and boundary layer. 


up to 80 keV (5 keV and 10 keV wide) were obtained us- 
ing the Earth occultation technique. The HR was defined as 
20-40 keV/ 10-20 keV and we used the flux in the 10-80 keV 
band for the intensity in our diagram. To reduce large uncer- 
tainties, the light curves were rebinned and then the HRs were 
computed. It is to be noted that like all HRs, these are instru- 
mentally dependent. From that figure we can see that the vast 
majority of observations reveal a constant HR. Only at low flux 
intensity levels there seems to be a smooth hardening of the 
spectrum of A 0535+26 as the flux intensity increases. There- 
fore, the diagonal and horizontal branches in the HRI diagram 
observed by Reig (2008) in a correlated spectral-timing analy- 
sis of four BeXRBs are not present for A 0535+26. We point 
out that the study by Reig (2008) was carried out during major 
X-ray outbursts and not all sources displayed the entire pattern 
of variability. 

3.3. Optical Measurements 

The middle and bottom panels of Figure 1 show the Ha 
EW and V magnitude long-term evolution for HDE 245770 
since ~ 1978. The compilation of Ha EW measurements come 
from Aab (1985, stars), Clark et al. (1998, filled circles), Coe 
et al. (2006, filled up triangles), Grundstrom et al. (2007, down 
triangles), and Moritani et al. (2010, squares). The historical 
observations of the V magnitude come from Gnedin et al. 
(1988, stars), Lyuty & Zaitseva (2000, up triangles), Zaitseva 
(2005, crosses), and Coe et al. (2006, right triangles), as well 
archived observations from HIPPARCOS (circles), ASAS3 
(plus), and INTEGRAL/ OMC (left triangles). Table 1 shows 
recent measurements of the Ha EW and the V magnitudes 
obtained with LT, NAOC, OSN, and OAO. 

In Figure 1 we see that, in the long term, both the Ha EW 
and the V magnitude show a parallel global trend. From MJD 
~44000 to 50000 both magnitudes presented an overall decreas- 
ing trend, and after the minimum is reached they climbed back 
up to 2009 November (~MJD 55136). Then, both magnitudes 
peaked in this period, and afterward a dramatic decrease took 
place, which continued through 2011 March. Furthermore, the 


shape of the Ha profile was more simple around 2009, evolving 
from a clear double structure (Coe et al. 2006) to a single peak. 
Later on, a double-peaked structure is observed (see Figure 4). 

4. TIMING ANALYSIS 

We have performed distinct types of timing analysis for 
A 0535+26 in X-rays. These include a study of the evolution 
of the pulse frequency, a search for QPOs, and a pulse profile 
study. 

4.1. Long-term Pulse Frequency Evolution 

The top panel of Figure 5 shows a long-term pulse frequency 
history for A 0535+26 obtained by several authors since 1975 
to 2006. Among them are Rosenberg et al. (1975), Fishman & 
Watts (1977), Bradt et al. (1976), Li et al. (1979), Hameury et al. 
(1983), Frontera et al. (1985), Nagase et al. (1982), Sembay 
et al. (1990), Motch et al. (1991), Coe et al. (1990), Makino 
(1989), Sunyaev & Titarchuk (1989), Cusumano et al. (1992), 
and Finger et al. (1994, 2006, 2009a). 

A 0535+26 shows a brief spin-down period (MJD 
42000-44000) followed by a global spin-up trend (MJD 
44000-49500), and ending with a spin-down trend. Although 
between MJD 50000 and 53500, no X-ray outburst activity was 
reported from this source (see Figure 1), Coe et al. (2006) de- 
tected a 110 ± 0.5 day orbital modulation pattern during this 
long-quiescent X-ray state with RXTE/ ASM. Hill et al. (2007), 
Negueruela et al. (2000), and Orlandini et al. (2004) detected 
pulsations at ~ 103.5 s with different instruments during this 
period. Then, suddenly the 2005 giant outburst took place and 
A 0535+26 exhibited a rapid spin-up. The following giant out- 
bursts in 2009, 2010, and 2011 allowed A 0535+26 to slowly 
recover from the overall spin-down period. We note that the 
transitions between giant outbursts resulted in global smooth 
spin-down periods. Nevertheless, in spite of this overall picture 
A 0535+26 might be just exhibiting a random walk in pulse 
period over the long term. 
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Wavelength (Angstroms) Wavelength (Angstroms) 

Figure 4. Recent Her line profile evolution from 2009 October up to 201 1 March. All spectra have been normalized by dividing by the mean flux calculated between 
6450 A and 6500 A. As a reference, some of the most recent outburst peaked around 2009 August 7, December 17 (giant), 2010 March 28, July 1, September 1, and 
201 1 February 28 (giant). 


4.2. Pulse Frequency Evolution Monitored by Fermi /GBM 

The middle panel of Figure 5 shows the recent pulse frequency 
history by GBM. We have used the orbital elements from 
Finger et al. (1996) with adjusted epoch and period, i.e., 
P orb = 111.07(7) days and 7 pl;nasll( , n = MJD 53613.00(16). After 
the first normal outburst detected by GBM in 2008 October, 
A 0535+26 experienced an overall steady spin-down trend up 
to the 2009 December outburst. Then, a rapid spin-up occurred 
followed by a return to a steady spin-down trend. This seems 
to be a common pattern for A 0535+26, also seen surrounding 
another giant outbursts. The bottom panel of the same figure 
shows the 12-25 keV pulsed flux for the same period by 
GBM. Once again, we find a strong pulsed flux/spin-up rate 
correlation, indicating the presence of an accretion disk around 
the neutron star. 

4.3. Quasi-periodic Oscillations in the 2009 
December Giant Outburst 

During the 1994 giant outburst, a broad QPO from 27 to 
72 mHz was detected with BATSE (Finger et al. 1996) in 


the 25-60 keV band. With GBM we analyzed the aperiodic 
variability in the X-ray flux of A 0535+26 following the 
method described in Finger et al. (1996). During the 2009 
December giant outburst we found that the power spectra of the 
hard-X-ray flux of this source, between 3 mHz and 1 Hz, 
consisted of an approximately 1// power-law continuum plus 
a broad QPO and a pulse component (Finger et al. 2009b). On 
December 10, the QPO was centered on 62 ± 1 mHz (FWHM 
of 29 ± 2 mHz) in the 25-50 keV band. Figure 7 shows a 
color-scale representation of the evolution of the power spectra 
over the course of the outburst in the 25-50 keV band. The QPO 
is evident during much of the outburst, and shows a marked 
evolution in centroid frequency. We were able to detect the 
QPO from December 4 to 27. In Figure 6 we can see that this 
QPO was stronger in the 50-100 keV band, but not detected in 
the 12-25 keV range. 1 ' In the bottom-left panel of Figure 7, it 
is also shown how the center frequency rose from 30 to 70 mHz 


17 We note that this is the most sensitive energy band in GBM. Further studies 
investigating these results will be carried out in a forthcoming paper using data 
from other X-ray missions, as for instance RXTE and Suzaku. 
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Figure 5. Top: long-term frequency history of A 0535+26 since 1975 (see the text for references). Middle: zoom of the frequency history of this source since 2008 by 
GBM. The frequencies are connected with a dot-dashed line. Bottom: daily average 12-50 keV pulsed flux measured with GBM. 

(A color version of this figure is available in the online journal.) 


and went back to 30 mHz. In addition, a strong QPO center 
frequency /X-ray flux correlation is found (bottom-right panel 
of Figure 7). When we perform a linear fit to the QPO center 
frequency and X-ray flux data sets, we obtained for the best fit 
a correlation coefficient of 0.9677762. 

4.4. Pulse Profile Overview 

During the bright phases of an outburst, the pulse profile of 
A 0535+26 displays a complex shape in the soft-X-ray range, 
below 8 keV, and a simple double-peaked profile, with very 
different spectral shape of the pulses, at higher energies — quite 
similar to what is observed in several other accreting pulsars, 
e.g., SAX 2103.5+4545 (Camera Arranz et al. 2007). A similar 
overall shape has been found in a wide range of observations 
over several decades in time and outbursts of quite different 
peak brightness (e.g., Bradt et al. 1976 and more recently 
Kretschmar et al. 2005, Finger et al. 2006, Camero-Arranz 2007, 
and Caballero et al. 2007). 

Using OSSE data, Maisack et al. (1996) obtained pulse 
profiles of A 0535+26 in four energy bands (from 35 to 100 keV) 
during the first and the second half of the 1994 giant outburst. 
They saw that the intensity increase between the first and the 
second half of the outburst was most pronounced at lower 
energies, indicating that the spectrum becomes softer as the 
overall intensity rose. The pulse profile shape was similar to 


that reported from the previous giant outburst in 1989 with 
HEXE (Kendziorra et al. 1994; Kretschmar et al. 1996). 

Bildsten et al. (1997) confirmed this behavior based on 
BATSE data of the 1994 giant outburst. Moreover, they found 
that at low X-ray luminosities (^1.6 x 10 36 erg s') the 
double-peak pulse profile structure evolved into a broad sin- 
gle sinusoidal-like peak, becoming a top-hat shape at higher 
energies. 

4.4.1. Recent Giant Outbursts 

Using GBM data we obtained for the first time a continuous 
coverage of the pulse profile evolution for A 0535+26 during a 
giant outburst at low energies, in the 8-12 keV band. Figure 8 
shows the pulse profile evolution with energy for an observation 
at the peak of the 2009 December giant outburst. The GBM 
data show a double-peaked pulse profile structure with a large 
dip between the pulses (two cycles are shown for clarity). The 
first component spans phases from 0 to ~0.3 and the second 
one from ~0.5 to 1, with the dip located between phases ~0.3 
and ~0.5. In that observation we can also see how the two main 
components of the profile evolve in opposite ways with energy. 
With increasing energy, the profile switches from a weak first 
and strong second peak to a strong asymmetric primary and 
almost vanishing secondary peak. In our study, we found that 
in the bright outbursts of 2010 and 2011 the strength of the two 
peaks is more balanced, albeit at a much lower peak luminosity. 
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Table 1 

Recent Her EW (Error ±0.5) and V Magnitudes 


MJD 

Ho- EW 

V mag 

Error 

Telescope 

55060.50 


8.900 

0.05 

OAO 

55082.00 

-20.7 



LT 

55083.00 

-21.3 



LT 

55084.50 

-21.6 



LT 

55087.00 

-21.5 



LT 

55087.50 

-21.0 



LT 

55091.00 

-21.0 



LT 

55115.50 


9.120 

0.01 

OSN 

55116.50 


9.110 

0.01 

OSN 

55117.50 


9.100 

0.01 

OSN 

55127.88 

-24.70 



NAOC 

55128.50 


9.150 

0.01 

OSN 

55128.84 

-24.82 



NAOC 

55129.50 


9.140 

0.01 

OSN 

55129.83 

-24.38 



NAOC 

55130.82 

-25.58 



NAOC 

55131.50 


9.150 

0.01 

OSN 

55132.85 

-25.02 



NAOC 

55133.82 

-25.86 



NAOC 

55145.00 

-23.6 



LT 

55147.00 

-23.5 



LT 

55150.00 

-23.6 



LT 

55154.00 

-24.6 



LT 

55158.00 

-24.7 



LT 

55162.00 

-24.0 



LT 

55166.00 

-24.5 



LT 

55168.00 

-23.9 



LT 

55171.00 

-23.4 



LT 

55175.00 

-23.3 



LT 

55176.81 

-23.14 



NAOC 

55191.00 

-23.6 



LT 

55196.00 

-23.7 



LT 

55199.00 

-24.2 



LT 

55203.00 

-23.8 



LT 

55208.00 

-23.5 



LT 

55211.00 

-23.2 



LT 

55213.00 

-21.9 



LT 

55214.50 


9.230 

0.01 

OSN 

55218.00 

-22.2 



LT 

55220.00 

-22.6 



LT 

55223.90 

-23.0 



LT 

55309.00 

-22.0 



LT 

55309.90 

-22.0 



LT 

55487.00 

-11.0 



OSN 

55494.00 

-11.8 



LT 

55495.10 

-11.8 



LT 

55496.10 

-12.1 



LT 

55499.20 

-11.9 



LT 

55502.10 

-10.4 



LT 

55505.10 

-10.3 



LT 

55508.20 

-10.7 



LT 

55513.10 

-11.5 



LT 

55516.00 

-12.1 



LT 

55517.10 

-12.3 



LT 

55520.20 

-12.4 



LT 

55523.10 

-11.1 



LT 

55542.42 


9.269 

0.02 

OSN 

55544.00 

-10.5 



LT 

55545.10 

-10.5 



LT 

55545.52 


9.285 

0.02 

OSN 

55546.44 


9.293 

0.02 

OSN 

55554.20 

-11.3 



LT 

55564.00 

-10.2 



LT 

55565.10 

-10.2 



LT 

55572.90 

-9.3 



LT 


Table 1 

(Continued) 


MJD 

HcrEW 

V mag 

Error 

Telescope 

55574.00 

-9.2 



LT 

55574.00 

-9.5 



OSN 

55574.90 

-9.2 



LT 

55578.90 

-10.1 



LT 

55580.00 

-10.0 



OSN 

55580.10 

-10.5 



LT 

55595.52 


9.271 

0.02 

OSN 

55599.46 


9.272 

0.02 

OSN 

55600.48 


9.266 

0.02 

OSN 

55601.45 


9.269 

0.02 

OSN 

55602.44 


9.276 

0.02 

OSN 

55603.51 


9.268 

0.02 

OSN 

55604.44 


9.263 

0.02 

OSN 

55605.43 


9.265 

0.02 

OSN 

55608.90 

-10.0 



LT 

55610.90 

-11.0 



LT 

55614.90 

-12.0 



LT 

55622.90 

-10.8 



LT 

55664.90 

-9.2 



LT 



0.01 0.10 
Frequency (Hz) 


Figure 6. QPO centered at 62 mHz from a GBM observation on 2009 December 
1 1 in three energy bands. We can clearly see that this QPO is stronger in the 
50-100 keV band but not detected in 12-25 keV range. Dotted lines denote 
pulse harmonics. The daily mean pulse profile was subtracted from the data 
before power spectra was made. 

During these latter two outbursts the energy dependence of the 
pulse profile is in line with that observed, e.g., in the 1994 giant 
outburst by BATSE (Bildsten et al. 1997), i.e., two peaks of 
approximately same strength at lower energies with the second 
decreasing with increasing energy (see Figure 8). 

In the top panel of Figure 9, we show a three-dimensional 
pulse profile evolution with time and intensity for the 2009 
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Figure 7. Top: dynamic power spectra for A 0535+26 showing the evolution of the QPO center frequency in the 25-50 keV band. The color scale gives the fractional 
power per logarithmic frequency interval. Bottom: QPO center frequency evolution vs. time (left) and vs. the pulsed flux in the 25-50 keV band (right). A strong QPO 
center frequency /X-ray flux correlation is found. Red points indicate the rise of the outburst and black the fall. The best linear fit gives a correlation coefficient of 
0.9677762 (dashed line). 

(A color version of this figure is available in the online journal.) 


December giant outburst in the 12-25 keV range. Darker colors 
denote lower flux intensities (see the electronic version for color 
plots). We see that the two main components become more 
smooth and balanced as the flux intensity decreases, finally 
merging into one single component (almost flat) at very low 
fluxes, as previously seen by BATSE and/or in normal outbursts. 

4.4.2. Recent Normal Outbursts 

Following the evolution of the pulse profile with luminosity, 
at energies above 8 keV we found with GBM a broad single 
sinusoidal peak at lower luminosities and a clearly double-peak 
profile at high luminosities (see Figure 9). With RXTE / PCA we 
found a similar evolution as in previous studies for this source 
by different authors. With increasing luminosity the low-energy 
profiles became more complex, exhibiting multiple components 
while the high-energy profiles evolved from double to single 
peaks. RXTE data of the 2005 September outburst showed in 
addition that the pulse profile changed shape around the energy 


of the fundamental cyclotron line at 45 keV (Camero-Arranz 
2007). Moreover, the pulse shape behavior during the pre- 
outburst peak in that outburst was observed to be different than 
during the main peak (Camero-Arranz 2007; Caballero et al. 
2007). This has not been observed in other occasions. 

5. DISCUSSION 

5.1. Torque / Luminosity Correlations 

BATSE observed BeXRB peak spin-up rates of 4.3 x 
10~ n Hz s" 1 (2S 1417-624), 3.8 x 10" 11 Hz s" 1 (GRO 
J1750-27), and 8 x 10“ 12 Hz s" 1 (4U 0115+634; Bildsten 
et al. 1997). The discovery outburst of EXO 2030+375 found 
it spinning up at a rate of 2.2 x 10~ n Hz s 1 (Reynolds et al. 
1996), and for SAX J2103. 5+4545 a peak spin-up rate of 1.4 x 
10~ 12 Hz s -1 has been reported (Camera Arranz et al. 2007). 
For other type of accreting binaries, e.g., the accreting millisec- 
ond pulsars (AMSPs), there are only a few sources of this type 
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Figure 8. Left: pulse profile correlation with energy during the 2009 December giant outburst near the maximum (from top to bottom: 8-12, 20-25, 25-50, and 
50-100 keV bands) by GBM. Two cycles are shown for clarity. Middle: pulse profiles near the maximum of the large event in 2010 March (from top to bottom: 8-12, 
20-25, and 25-50 keV bands). Right: pulse profiles of the giant outburst in 201 1 March in the same bands. 


where a spin-up could be measured at a rate between 1 x 10” 13 
and 1 x 10” 12 Hz s” 1 (see, e.g., Papitto et al. 2008), compatible 
with the expected accretion-driven spin-up. On the other hand, 
timing noise dominates the pulse phases observed from other 
AMSPs, among them SAX J 1808.4-3658, so that a measure of 
the spin evolution during outburst is an extremely difficult task 
(Patruno et al. 2012; Papitto et al. 2011). It is to be noted that the 
long-term spin-down for this source has remained stable over 
13 years with u ~ — 1 x 10” 15 Hz s” 1 , compatible with the 
magnetic dipole spin-down of a neutron star with a field of few 
10 8 G (Patruno et al. 2012). 

Previous studies on A 0535+26 showed the existence of 
a strong torque/X-ray luminosity correlation in this BeXRB 
system (see, e.g., Bildsten et al. 1997). Our results using GBM 
data are in very good agreement with those already published. 
In Figure 5 we can see a clear torque-flux correlation, which 
suggests disk accretion at least in the larger outbursts. We 
caution that the apparent frequency rise during the normal 
outbursts may be an artifact due to the uncertainty of the 
periastron epoch. The peak spin-up rate for the 1994 giant 
outburst was 11.9 x 10 12 Hz s 1 (Finger et al. 1996), 
comparable to the 2009 December giant outburst with a spin- 
up rate of 10.15(6) x 10” 12 Hz s” 1 . For the 2011 February 
giant outburst the peak of the spin-up rate was 6.36(5) x 


10” 12 Hz s” 1 . This value is smaller than in the previous giant 
outbursts, although we note that this event was weaker. 

Before BATSE, it was already known that the spin frequency 
in some transients decreases between outbursts (Bildsten et al. 
1997, and references therein). This was explained due to the 
propeller effect, when M becomes small enough so that the 
magnetospheric radius exceeds the corotation radius. Accretion 
is then centrifugally inhibited, and material may become tempo- 
rally attached to magnetic field lines, removing angular momen- 
tum and causing the star to spin down (Illarionov & Sunyaev 
1975). Finger et al. (1994) reported spin-down between out- 
bursts in A0535+26 at a rate of —2.2(6) x 10” 13 Hz s” 1 . More 
recently, before the 2009 December giant outburst, the spin- 
down rate between outbursts ranged from — 1 x 10” 13 Hz s” 1 to 
— 1.5 x 10” 13 Hz s” 1 . After this event, the spin-down between 
the last normal outburst and the last giant outburst occurred 
in 2011 February happened at an approximate rate of —3.3 x 
10” 13 Hz s” 1 . 

A 0535+26 also exhibits global spin-down trends during long 
periods of quiescence. Assuming that in quiescence A 0535+26 
enters the subsonic propeller regimen, it is expected to spin down 
at a rate of v = -4 jt u 2 /r 2 (GM)” 1 /” 1 ~ -2.4 x lO-^Hzs” 1 , 
where /i is the neutron star magnetic moment, M is the mass 
(1.4 M 0 ), and /is the moment of inertia (Henrichs 1993). Based 
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2 ° r 



Figure 9. Top: three-dimensional pulse profile evolution with time and flux intensity during the 2009 December giant outburst in the 12-25 keV band. Two phase 
cycles are shown. Darker colors denote lower flux intensities. The two main components of the pulse profile evolve into a single component (almost flat) at very 
low fluxes, as previously seen by BATSE. Bottom: three-dimensional pulse profile evolution with time and intensity during the double-peak normal outburst in 2009 
November, in the same band. 

(A color version of this figure is available in the online journal.) 


on our observations, during the quiescence period between 
2005 and 2008, A 0535+26 showed a spin-down trend with 
v ~ — 1.01(5) x 10- 13 Hz s _1 . These two values are not in 
agreement and therefore this behavior cannot be well explained. 
On the other hand, mass prevented from accreting by the 
presence of a centrifugal barrier does not necessarily have to 
leave the system (D’Angelo & Spruit 2012, and references 
therein). These authors found that the accretion from a disk 
on to a magnetosphere tends to take place in “trapped” states, 
where the accretion can be steady or cyclic. Two forms of cyclic 
accretion are found in this study, (1) the accretion-dead disk 
cycle (Syunyaev & Shakura 1977) and (2) where a pile-up of 
mass at the centrifugal barrier outside corotation is followed by 
an episode of accretion emptying the pile and the beginning of a 
new cycle. In this model, the effect of these cycles is to increase 
the average spin-down torque on the star. However, whether or 
not this model may be applicable to A 0535+26 has not been 
explored. 

5.2. Quasi-periodic Oscillations 

Reig (2008) carried out a rapid spectral and timing variability 
study of BeXRBs during type II outbursts. In this study, the 
aperiodic variability of the four BeXRBs was described with a 


relatively low number of broad Lorentzian components, and 
with peaked noise (for some of the sources as QPO noise) 
between the low-frequency and the high-frequency noise. In 
addition, the flat-topped noise found at lower frequencies for 
sources in the horizontal branch of the HRI diagram turned 
into power-law noise in the diagonal branch. Some of these 
components can be associated with the same type of noise as 
seen in low-mass X-ray binaries, although the characteristic 
frequencies are about one order of magnitude lower (Reig 2008, 
and references therein). 

The origin and properties of QPOs are poorly understood. In 
accretion-powered X-ray pulsars, QPOs fall in the mHz range. 
The frequencies of these QPOs range from 1 to 400 mHz, a 
few orders of magnitude lower than in low-mass X-ray binaries 
(Reig 2008, and references therein). It is generally accepted 
that QPO features correlated with flux/spin-up are a powerful 
indication of the presence of an accretion disk. According 
to the beat frequency (BF) model (Alpar & Shaham 1985; 
Lamb et al. 1985), blobs of matter orbiting at the inner disk 
boundary at approximately the Keplerian frequency (KF) v k 
are gradually removed through the interaction with the neutron 
star magnetosphere. The magnetic field rotates with the neutron 
star at a frequency v s , modulating the accretion rate and source 
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luminosity. Therefore, the QPO frequency results from the beat 
between the orbital frequency of the blobs at the inner disk 
boundary and the spin frequency v s of the neutron star (uqpo = 
Vk - i> v ). In the KF model (van der Klis et al. 1987), the inner edge 
of the accretion disk contains structures that persist for a few 
cycles around the neutron star, and modulate the observed flux 
by obscuration producing a QPO at the Keplerian frequency i+. 
At lower mass accretion rates, torques caused by the interaction 
of the magnetic field with the accretion disk also need to be 
considered. In the magnetically threaded disk model (Ghosh & 
Lamb 1979), a magnetic torque is generated by the neutron star 
magnetic field threading the disk. In their model, the coupling 
between the neutron star and the accretion disk occurs in a broad 
transition zone located between the flow inside the neutron star 
magnetosphere and the unperturbed disk flow. Another model, 
the magnetic disk precession model (Shirakawa & Lai 2002), 
attributes the mHz QPO in X-ray pulsars to warping/precession 
modes induced by magnetic torques near the inner edge of the 
accretion disk. We note that this model has only been applied to 
weakly magnetized neutron stars in low-mass X-ray binaries. 

For A 0535+26, Finger et al. (1996) detected for the first 
time an mHz QPO during the 1994 giant outburst observed 
by BATSE. This QPO was detected for 33 days during the 
outburst. During this time the center frequency of the QPO 
rose from 27 to 72 mHz and then fell again to 25 mHz. These 
observations showed a strong correlation between the X-ray 
flux and spin-up rate for A 0535+26, reinforcing the idea of 
the presence of an accretion disk (Finger et al. 1996, and 
references therein). These authors successfully applied both 
the BF and the Keplerian models to these observations. These 
results showed that the QPO center frequency was, like the 
spin-up rate and the flux, controlled principally by the mass 
accretion rate and linked with the inner region of the accretion 
disk. In addition, Bozzo et al. (2009) discussed the physical 
equations that determine the position of the inner disk boundary 
by using different prescriptions for the neutron-star-accretion- 
disk interaction. They applied the results to several accreting 
pulsars showing both QPOs in their X-ray flux and torque 
reversals. They also concluded that for A 0535+262 a good 
agreement between the BF model and the Ghosh & Lamb (1979) 
model was obtained. 

In the present work, we reported on the detection of a hard 
30-70 mHz QPO using GBM data for A 0535+26, during the 
giant outburst occurred in 2009 December (Finger et al. 2009b). 
The QPO was detected for 24 days and the center frequency 
rose from 30 to 70 mHz and went back to 30 mHz. Once 
again, as in the QPO detected by BATSE, a strong QPO center 
frequency /X-ray flux correlation was found. The flux QPO 
frequency correlation may be explained with the BF model. 
However, the fact that the QPO was not detected by GBM in 
the most sensitive band of 12-25 keV, but only above this one, 
is inconsistent with this interpretation. It is outside the scope 
of this work to elaborate on a sophisticated theory reproducing 
these observations, however we will speculate with the origin 
of this inconsistency. In the KF model, the QPO was explained 
as produced by QPOs in the thick inner accretion disk, which 
partially obscures the radiation emitted in our line of sight. We 
can identify the major contribution of the hard-X-ray component 
in the spectrum of A 0535+26 with emission from the lower 
part of the accretion column and from the central disk regions, 
and the low-energy component with emission from scattering 
from the neutron star surface. One possible explanation for our 
observations is that the regions where hard X-rays originate 


might be obscured in our direction, while for the softer X-ray 
radiation no QPO is produced because any appreciable emitting 
region is occulted. 

QPO studies may also allow us to constrain the accretion 
flow. Considering the BF model and the relationship between 
the vqpo and the Keplerian radius ( r k = (GM/4 k 2 Vj) 1 / 3 ), we 
obtained a value for r ^ of ~9100 km (~9960 km using the KF 
model), for an observation at the peak of the 2009 December 
giant outburst. For A 0535+26 the last detection of a QPO by 
GBM was in 2009 December 27, with a vqpo of ~33.771 mHz. 
This means an r k ~ 1 3,500 km (with KF model ~ 16,060 km), 
consistent with the standard picture of an accretion disk moving 
outward as the accretion rate decreases. 

Finger et al. (1996) concluded that the QPO in A 0535+26 
is consistent with the BF model, with the QPO occurring at the 
beat frequency vqpo = i+ — v ns . They however believed that 
the fundamental cyclotron line energy was near 110 keV. This 
is now known to be the first harmonic, with the fundamental 
line near 50 keV (Kendziorra et al. 1994). We therefore will 
reconsider this interpretation. 

This model predicts two relationships between observed 
quantities (Finger et al. 1996). The first is between the pulsar 
spin-up rate v ns and the Keplerian frequency v/ f , 

v ns = X N v 2 with Bat = 1 .6K 1/2 B 2 R 6 (GM)~ l I~\ (1) 

where B is the surface magnetic field at the magnetic pole, R 
is the neutron star radius, G is Newton’s gravitational constant, 
M is the neutron star mass, and I is its moment of inertia. K 
is the ratio of the inner disk radius ro to /r 4 / 7 (GM) -1 / 7 M~ 2 ^ 7 , 
which is expected to be near 1 (in particular, K ~ 0.91 gives 
the Alfven radius for spherical accretion). Here, /i — 0.5 BR 3 
is the magnetic moment, and M the mass accretion rate. The 
value determined from the 1994 giant outburst was B v = 
(1.91 ±0.02) x 10“ 9 ergcm- 2 s 4/3 . 

The second relationship is between the flux F and the 
Keplerian frequency v*, 

F = X F vl /3 with B f = lAafid- 2 K 1/2 B 2 R 5 (GM)~ 2/3 , (2) 

where a is the fraction of the bolometric luminosity in the 
observed energy band, accounts for beaming, and d is the 
source distance. From the peak QPO frequency of 69 ± 1 mHz 
in the 2009 giant outburst, and the 3-50 keV peak flux 
of 1.25 x 10~ 7 erg cm -2 s _1 , we find B/r = (4.7 ± 0.2) x 
10~ 5 ergcirT 2 s 4/3 . 

If we determine B — 4.3 x 10 12 G from the cyclotron line 
energy of 50 keV, and adopt the nominal values a = ft — l,d = 
2 kpc, R = 10 6 cm, and M — 1.4 M 0 , then from Equation (2) 
and the measured value of B/r we would conclude that K — 2.4. 
This would imply that the QPO is not formed at the inner region 
of the accretion disk but at a radius much further out, or that the 
flow at the inner region of the disk is sub-Keplerian. We point 
out that the resulted discrepancy of the value of K between 
the 1994 ( K ~ 1) and 2009 giant outbursts is due to different 
assumptions, in particular, the implied magnetic field and the 
assumed values of a, /L and d. In Finger et al. (1996), they 
proposed values fora ss 0.5, 2 ^ j3 ^ 3, 1.3 kpc ^ d ^ 2.6 kpc 
and, as stated before, with the fundamental of the cyclotron line 
energy at 1 10 keV. 

The other possibility is that the adopted values for the 
parameters needed to convert the measured flux to the true 
luminosity, and the neutron star parameters, can be adjusted in a 
plausible way so that the BF model agrees with observations. In 
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Figure 10. Combination of neutron star parameters in Equations ( 1 ) and (2) are 
plotted for the set of neutron star models tabulated in Wiringa et al. (1988). The 
TI curve is for the tensor interaction model of Pandharipande & Smith (1975), 
which is a very stiff equation of state. The PA curve is for the Pandharipande 
(1971) hyperon model, which is a very soft equation of state. The AV14+UVII, 
UV14+UVII, and UV14+TNI models are from Wiringa et al. (1988). The 
symbols on the curves give the points where the mass is 0.5 A/q (triangles), 
1.0 AJq (squares), and 1.4 Mq (diamonds). The vertical dotted lines are for 
constant values of K, which determines the ratio of the inner disk radius to the 
Alfven radius (K 0.91 gives the Alfven radius for spherical accretion). The 
sloped dotted lines are for constant values of d(af})~ l t 2 , where d is the source 
distance, a is the bolometric correction factor, and /I is the beaming factor. 

Figure 10 we show the combinations of neutron star parameters 
in Equations (1) and (2) for a set of neutron star models. The 
plot also shows lines of constant K and d(apy' /2 assuming 
a magnetic field of B = 4.3 x 10 12 G with the measured 
values of and R F . The conversion of the measured flux 
to the true luminosity depends on d(a/i i) -1 / 2 . The most recent 
measurement of the source distance d is ~2 kpc (Steele et al. 
1998), with previous measurements ranging from 1.3 to 2.6 kpc. 
The bolometric fraction a probably above 80%. The beaming 
factor p result from the fact that the phase averaged flux only 
samples the pulsar’s spherical beam pattern over a cone. This 
term may vary from 1 with viewing direction by ~50%. The 
factor K is expected to be in the range of about 0.5-1 .0. So overall 
we think the BF model is plausible for neutron star parameters 
in the region with K in the range 0.5-1, and d(afi)~ l G in the 
range 1-3 kpc. This contains neutron stars with stiff equations 
of state, and masses of 1 M 0 or less. We note, however, that 
the current mass estimates for neutron stars are in excess of 
1.2 M q , with the majority of normal pulsars lying in the mass 
range 1.4-1. 8 M 0 (Ozel et al. 2012; Steiner et al. 2010). 

Reynolds & Miller (2010) found broadened Fe K„ profiles in 
two Chandra grating observations that occurred during the tail 
of the 2009 giant outbursts. They fit the lines with Gaussians 
with sigmas of ~ 5 1 00 and 4100 km s _1 , respectively. They 
interpreted broadening as due to the azimuthal velocity at the 
inner edge of the accretion disk. These observations occurred 
at approximately 50% and 40% of the peak flux, respectively. 
If the BF model for the QPO is correct, we can predict the KF 


at the disk inner edge, which scales as the F 3 / 7 , giving v* = 
60 and 53 mHz, respectively, for the two Chandra observations. 
The azimuthal component of velocity at the inner edge is then 
or 3700 and 3500 km s -1 . Given that 
these set the maximum amplitude of the projected velocity 
distribution, and these will be decrease further by the (unknown) 
inclination of the disk from the line of sight, we find this 
prediction in disagreement with the observations. 

5.3. Pulse Profile Study 

Our study is mostly in line with previous observations (see 
Section 4.4). At medium to high source luminosity, we find 
a multi-peak pulse shape at energies below ~8 keV, evolving 
into a double-peak structure at higher energies. The high-energy 
pulse shape evolution with luminosity also agrees with previous 
works, in that the profiles evolve toward a single broad peak as 
luminosity decreases (see Figure 8). 

The giant outburst of 2009, the brightest in our sample of 
outbursts monitored by GBM, is peculiar in that the two main 
components of the pulse profile evolved in opposite ways with 
energy. To our knowledge, this has been never reported before. 
In contrast, for the other bright outbursts in 2010 and 2011, 
our results are in good agreement with the energy dependence 
seen in the 1994 giant outburst (Bildsten et al. 1997). At the 
moment we have no explanation for the differences between 
the outbursts. The mechanisms leading to different kinds of 
outbursts are not well understood for different sources and 
currently there is no good model to describe A 0535+26. It 
remains unknown whether the giant outburst detected by BATSE 
in 1994 was in fact different from the giant outburst detected 
by GBM in 2009. With BATSE the pulse profile study was 
carried out above 20 keV, preventing us from performing any 
comparison. Above 20 keV the evolution of the pulse profile is 
consistent with GBM. 

Interpreting pulse profiles in terms of the underlying emission 
pattern is far from evident since various effects, including highly 
anisotropic cross-sections for emission and scattering and light- 
bending around the neutron star must be taken into account 
(Kraus et al. 2003). For A 0535+26, Caballero et al. (201 la) have 
disentangled the contribution of the two single pole components 
for A 0535+26 using pulse profiles from the 2005 August 
normal outburst. The beam pattern was interpreted in terms 
of a geometrical model that consists of a hollow column plus 
a halo of scattered radiation around the neutron star surface. In 
the model by Ibragimov & Poutanen (2009), the presence (or 
not) of a secondary maximum in the pulse profile is interpreted 
by the appearance (or disappearance) of a secondary antipodal 
spot. This model explains that obscuration of the antipodal spot 
by variations of the inner disk radius reproduces the variety of 
the pulse profiles observed in the accreting millisecond pulsar 
SAX J1808. 4-3658. The applicability to A 0535+26 is not yet 
certain. 

An intriguing case was the pulse shape behavior at the rise 
of the 2005 September normal outburst (Camero-Arranz 2007; 
Caballero et al. 2007). One of the components of the multi- 
peak structure at low energies became more prominent com- 
pared with the pulse profile during the main outburst. The 
evolution with energy was also different, evolving into one 
strong simple peak as the energy increased. A different accre- 
tion regime was suggested during the pre-outburst phase than 
during the main outburst. A possible explanation was bursts of 
wind-like accretion from capture of low angular momentum 
material, while the higher angular momentum material forms 
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Time (MJD ) 

Figure 11 . Zoom of Figure 1 from ~2002 August up to 201 1 March. We can see the anti-correlation of the Ha EW and the visual magnitude, with each giant X-ray 
outburst occurring during a decline phase of the optical brightness. 

(A color version of this figure is available in the online journal.) 


an accretion disk, which feed the main part of the outburst. 
Caballero et al. (2007, 2008) found a change of the cyclotron 
line energy up to ~50 keV during this period and no significant 
change in the spin frequency up to the time the main outburst 
began. Postnov et al. (2008) proposed that magnetospheric insta- 
bilities developing at the onset of the accretion were the origin 
of the flaring activity. This may explain the difference in the 
cyclotron line energy and energy-dependent pulse profiles that 
remained stable as the neutron star starts spinning up and the 
flaring activity disappears. We point out that for A 0535+26 
while the profile changes in the later stages of an outburst the 
line actually seems stable (Caballero et al. 2010). 

In contrast, during the double-peaked outbursts of 2009 
August and 2010 July (see Figure 2) no changes in the shape of 
the pulse profile were found in the present work. Also no change 
of the cyclotron line energy was found by Caballero et al. (2009) 
for the former outburst. We leave as an open question whether 
these pre-outbursts peaks are related to the formation and/or 
evolution of the disk around the Be star. 

5.4. Be-disk-Neutron-star Interaction 

In this study we have shown how the Ha line is practically 
always in emission, except for a small period around MJD 


51000, telling us about the presence of a persistent but variable 
Be disk. In the long term, the optical brightening of the Be star, 
besides large negative values of the Ha EW, seems to precede 
the apparition of episodes of high X-ray luminosity (and neutron 
star spin-up), including giant X-ray outbursts, indicating that 
the variations in the circumstellar structure are the precursors 
of this X-ray activity. This behavior has been observed in other 
BeXRBs, as for instance in EXO 2030+375 (Wilson et al. 2008). 
We have also observed that the onset of both normal and giant 
outbursts are initiated near periastron. Although this does not fit 
the Stella et al. (1986) picture of type II outbursts with timing 
unrelated to the underlying orbital period. 

Taking a closer look, the Be star lost and renewed its 
circumstellar disk during the period from ~1996 until 2005. 
After this period the Ha EW and the optical brightening were 
no longer correlated. Figure 1 1 shows that after ~MJD 52500 
(2002 August) each giant X-ray outburst occurred during a 
decline phase of the optical brightness. It is interesting that the 
Ha showed a strong emission, while its optical brightness had a 
decreasing trend before and after the 2009 giant X-ray outburst. 
That is, the Ha EW and the V-band brightness showed an anti- 
correlation. Yan et al. (2012) propose that this anti-correlation 
may indicate that a mass ejection event took place before the 
2009 giant X-ray outburst. 
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The Ha EW measurements (see, e.g., Haigh et al. 2004; 
Grundstrom & Gies 2006) implies this was the largest Be disk 
observed in the system since its discovery. This may explain the 
unusually active period of A 0535+26 in X-rays. The decreasing 
trend and drastic optical spectroscopic variability after this 2009 
December giant outburst has also been recently reported by 
Moritani et al. (2011). However, their measurements of the 
Ha EW before and during the 2009 giant outburst are much 
lower than our values and those from Yan et al. (2012). Most 
likely these differences are not real, but the product of two 
different instrument systems (spectral resolution and continuum 
determinations). 

5.4.1. Recent Ha Line Profile Evolution 

It is interesting to compare (Figure 4) the Ha line profiles 
of the disks before, during, and after the giant outbursts in 
2009 December (left hand panel) and 2011 February (right 
hand panel), as observed by the FRODOSpec spectrograph 
( R ~ 5000) on the LT. These data present a far higher time 
resolution (~ few days) than that typically obtained around the 
time of previous bursts (~1 month). 

Throughout the 2009 outburst the Ha profile was consistently 
single peaked with a slightly asymmetric profile characterized 
by a broader blue wing, and no significant change can be seen 
around the time of outburst. A comparison with the higher 
resolution ( R ~ 30,000-60,000) but more infrequent spectra 
presented by Moritani et al. (201 1), over the same time period, 
shows the same basic structure, although they are able to resolve 
the peak which shows a double peak with velocity separation 
of ~ 100 km s . On occasions small additional absorption and 
emission features are also seen in their spectra. We note that 
Naik et al. (2012) did not find any major changes in a series 
of JHK spectra taken during the X-ray quiescent phase and the 
giant X-ray outburst in 201 1 February-March. 

In late 2010 (prior to the 2011 February outburst) the disk 
can be seen to be in a very different state, with an Ha EW of 
only half that seen in late 2009 and early 2010. The profile is 
double peaked with a separation corresponding to ~270 km s _1 . 
Most striking are the rapid V/R (violet/red) variations in the 
strength of the two peaks. These have a period of ~25 days and 
can be seen to repeat over two cycle from late 2010 October 
until mid 2010 December. This is far shorter than the typical 
~1000 day periods seen in isolated Be stars (Hummel 1998). 
V/R variations with a timescale of ~1 yr were noted from 
A 0535+26 in the period following the giant outburst of 1994 
February/March by Clark et al. (1998; from 1994 
September until late 2006), who attributed them to the pres- 
ence of a global one-armed oscillation (Okazaki 1991). From 
mid 2010 December the profile changes again, with the blue 
peak remaining dominant and the red wing broadening as 
a red component moves to higher velocities. No specific 
change at the time of the 2011 February outburst is obvi- 
ous, with the red wing simply continuing its move to higher 
velocities. 

A general pattern of behavior may start to be inferred 
from the behavior of the system following the 1994 March 
and 2009 December giant outbursts, where the disk becomes 
weaker and shows V/R variability commencing ~6 months 
following a giant outburst. Following the model of Okazaki & 
Negueruela (2001; which is consistent with the variability of 
the circumstellar disk after the 2009 giant X-ray outburst), it is 
important to note that those disk changes are not caused by the 
outburst, but rather that the outburst is a signal that changes in the 


disk are underway, presumably due to changes in the truncation 
resonance radius. If this pattern is to be believed, then when 
the source becomes observable again in 2011 September, we 
would expect the disk strength and associated line profile to 
differ from that presented in 201 1 April. In this overall picture 
what is still an open question is the mechanism/s that trigger a 
giant outburst. 

6. CONCLUSIONS 

A 0535+26 exhibits an unpredictable X-ray behavior, with 
long periods of quiescence related to an absence of Ha in 
emission, with periodic normal outbursts generally taking place 
before and/or after a giant outburst. Some peculiar normal 
outbursts of A 0535+26 present a double-peaked structure, with 
the origin of these pre-outburst peaks still unknown. Our study 
of the spectral variability of this source using GBM data from 
recent outbursts shows that there is a hardening of the spectrum 
as the flux intensity increases. 

Our pulse profile study for A 0535+26 confirms in general 
previous results. However, for the brightest outburst monitored 
by GBM in 2009 December, we found that the two main 
components of the pulse profile evolved in opposite ways 
with energy, which to our knowledge has been never reported 
before. A hard 30-70 mHz X-ray QPO was detected with 
GBM during this giant outburst. The new insights are that 
it becomes stronger with increasing energy and disappears at 
energies below 25 keV. Current models cannot explain these 
observations. 

In spite of the fact that in the long term a strong 
optical/X-ray correlation has been found for this BeXRB, how- 
ever, it is still difficult to foresee when a giant outburst is go- 
ing to take place and therefore what the trigger mechanism 
might be. On the other hand, in the medium term the Ha EW 
and the V-band brightness show an anti-correlation after ~2002 
August (MJD 52500). Each giant X-ray outburst occurred dur- 
ing a decline phase of the optical brightness, while the Ha 
showed a strong emission. We would like to point out that the 
large event that came after the 2009 December giant outburst 
does not follow the standard classification as either type I or 
type II outburst. 

The Be circumstellar disk was exceptionally large just before 
the 2009 December giant outburst. Most likely this was the 
origin of the unusual recent X-ray activity of this source. In 
late 2010 and prior to the 2011 February outburst, rapid V/R 
variations are observed in the strength of the two peaks of the 
Ha line, with a period of ~25 days, suggesting the presence of 
a global one-armed oscillation. A general pattern of behavior 
is proposed in this work, where the disk becomes weaker and 
shows V/R variability beginning ~6 months following a giant 
outburst. 
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